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The  present  computer  simulation  study  of  the  dislocation  distri¬ 
bution  near  the  surface  is  based  upon  a  modified  one-dimensional 
dislocation  model  which  considers  the  line  tension  of  a  disloca¬ 
tion  loop,  and  the  image  forces.  If  the  dislocations  are  gener¬ 
ated  from  a  surface  source  and  move  by  thermal  activation  over 
randomly  distributed  short  range  barriers,  an  inverse  pile-up  is 
created.  A  dislocation  density  gradient  exists  up  to  a  depth  of  __p> 
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C-jiO  pm)  for  a  total  of  tw$i)ty  dislocations  generated.  During  un¬ 
loading,  the  last  few  dislocations  generated  leave  the  crystal 
rapidjly,  and  within  three  hours  to  two  weeks  the  dislocation  den¬ 
sity  gradient  essentially  disappears  (depending  upon  the  para¬ 
meters  used).  This  is  in  agreement  with  the  experimental  results. 
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INTRODUCTION 


Whether  plastic  deformation  induces  a  soft  surface  or  a 
hard  surface  has  been  a  controversy  for  the  past  fifteen  years 
or  more.  The  experimental  evidence  of  a  higher  dislocation 
density  at  surface  or  in  the  core  region  following  plastic 
deformation  seems  often  contradictory.  There  has  been  sev¬ 
eral  mechanisms  proposed  which  would  lead  to  a  lower  dislocation 
density  in  the  near  surface  region  of  a  metal  with  a  clean 
surface.  The  object  of  the  present  study  is  to  develop  a  model 
which  will  result  in  a  higher  dislocation  density  in  the  near 
surface  region. 

The  present  model  is  based  on  two  well-known  facts : 

(1)  Surface  sources  can  be  activated  at  a  lower  stress  than 
sources  in  the  interior  of  a  sample.  (2)  Thermally  activa¬ 
ted  dislocation  motion  from  a  source  will  create  an  inverse 
pile-up.  Therefore,  if  all  dislocations  are  generated  from 
a  surface  source  and  their  motion  is  controlled  by  thermal 
activation,  then  there  will  always  be  a  higher  dislocation 
density  in  the  near  surface  region  during  loading. 

During  the  unloading  stage,  however,  the  dislocations 
will  travel  out  of  the  crystal  by  a  superposition  of  the 
following  factors:  (1)  Line  tension  due  to  the  curvature 
of  dislocation  loop  (2)  Attractive  force  by  image  disloca¬ 
tion  arrays,  which  are  mirror  images  to  the  real  dislocations. 
(3)  Repulsive  forces  by  those  dislocations  of  the  same  sign, 
lying  ahead  in  the  pile-up.  Therefore,  some  extra  barriers 


may  have  to  be  introduced  which  will  stop  the  dislocations 
from  exiting  too  easily.  The  intersection  of  a  secondary 
dislocation  slip  system  with  the  primary  ones  will  create 
these  extra  barriers. 

DISCUSSION  OF  RESULTS 

The  present  study  of  surface  dislocation  sources  will 
be  concentrated  on  the  following  three  models:  1.  A  one¬ 
dimensional  model  with  one  slip  system:  A  group  of  straight¬ 
edge  dislocations  parallel  to  each  other  are  created  from  a 
surface  source  and  move  on  a  slip  plane  by  an  applied  shear 
stress.  In  addition,  line  tension  terms  are  included  which 
take  into  account  the  circular  nature  of  dislocation  loops. 
Also,  the  image  forces  are  considered.  2.  A  two  dimensional 
model  with  two  slip  systems:  Two  systems  of  parallel- straight¬ 
edge  dislocations  lying  on  two  slip  planes  intersecting  at 
70°  which  can  operate  simultaneously  from  each  surface  source. 
3.  A  three  dimensional  model  with  two  slip  systems:  The 
static  interaction  between  various  configurations  of  two 
systems  of  nonparallel,  noncoplanar  straight  edge  disloca¬ 
tions  of  slip  system  [lrfl]  (l2l)  and  [Oil]  (2ll)  are  being 
considered. 

Isotropic  linear  elasticity  is  always  assumed  in  the 
present  study.  Only  room  temperature  is  considered. 
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Model  1 

The  thermally  activated  motion  of  straight  dislocations 
from  a  source  under  an  applied  stress  will  create  an  inverse 
dislocation  pile-up  on  a  slip  plane  with  randomly 
distributed  short  range  barriers  (SRB)  (Fig.  1).  The  equili¬ 
brium  position  of  a  dislocation  is  determined  by  a  static 
force  balance.  Normally,  under  an  applied  stress,  a  dislq  * 
cation  will  be  pushed  against  a  (SRB)  with  a  resultant  posi¬ 
tive  effective  stress  such  as  to  lower  the  activation  energy 
required  for  jumping  forward.  However,  there  are  occasions,  al 
though  very  few,  that  a  dislocation  "floats"  in  between  two 
SRB's  with  zero  effective  stress  and  thus  not  in  a  position 
to  jump  any  SRB  yet. 

In  order  to  relax  the  present  restriction  of  a  straight 
dislocation,  the  backward  force  component  of  the  line  tension 
of  a  dislocation  loop  is  considered.  If  the  source  is  at 
the  surface  of  the  crystal,  as  in  Fig.  2,  a  further  term  is 
needed,  i.e.,  the  algebraic  sum  of  the  image  forces  by  an 
array  of  inverse  pile-up  image  dislocations.  The  result  of 
adding  these  two  extra  backward  force  terms  adds  to  the  dif¬ 
ficulty  of  operating  a  dislocation  source. 

The  actual  distribution  of  dislocations  as  a  function 
of  time  depends  to  a  large  extent  on  the  chosen  values  of 
the  parameters  such  as  tq,  the  stress  for  dislocation  motion 
at  0  K,  t8,  the  source  operating  stress,  AGQ,  the  activation 
energy  at  zero  effective  stress,  and  finally,  AXav,  the 
average  spacing  among  randomly  distributed  SRB's.  In  the 


present  study,  because  the  exact  values  for  copper  are  not 
known,  some  reasonable  estimates  are  made  in  choosing  the 
appropriate  values  of  x^,  Tg  and  The  applied  stress, 

t  ,  is  chosen  so  as  to  give  a  reasonable  generation  time  for 
twenty  dislocations.  A  range  of  estimated  AGQ  values  were 
examined . 

It  is  found  that,  in  spite  of  the  differences  in  AGQ, 
the  penetration  depth  of  the  first  dislocation  generated  in 
an  array  of  twenty  is  always  between  500  and  600  pm.  This 
indicates  that  the  dislocation  distribution  is  approximately 
constant  for  a  given  number  of  generated  dislocations.  If  the 
dislocation  density  distribution  is  considered,  (Fig.  3)  the  density 
changes  rapidly  at  about  100  pm.  For  depths  greater  than 
100  pm  the  density  is  much  less  than  for  a  depth  up  to  100  pm. 

(This  is  in  agreement  with  the  experimental  results  of  Kramer 
who  has  stated  that  the  hard  layer  extends  into  a  depth  of 
approximately  100  pm.) 

As  soon  as  the  three  dislocations  nearest  to  the  surface 
leave  the  crystal  during  recovery,  the  dislocation  density 
established  by  the  inverse  pile-up  in  front  of  the  source 
essentially  disappears.  After  three  more  dislocations  leave 
the  crystal,  a  reversed  dislocation  density  gradient  is 
developed,  i.e.  a  lower  density  exists  near  the  surface. 

The  time  necessary  for  this  reversed  gradient  to  demonstrate 
itself  will  be  termed  "recovery  time".  The  dislocation  dis¬ 
tributions  after  recovery  are  similar  for  the  small  range  of 
AGq  values  compared  here.  (The  values  of  AGQ  choosen  represent 
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estimates  of  AGQ  for  copper) . 

As  can  be  seen  from  the  effects  of  different  AG_  in 

o 

Table  1,  both  the  generation  and  the  recovery  times  are  very 
sensitive  to  a  change  in  AGQ.  Therefore,  a  very  accurate 
experimental  determination  of  AGQ  for  copper  is  needed. 

The  present  study,  therefore,  points  out  the  possibili¬ 
ty  of  a  new  interpretation  of  the  controversy  between  "hard 
surface”  and  "hard  core"  theories  in  plastic  deformation. 

The  initiation  of  plastic  deformation  occurs  at  the  surface 
and  leads  to  higher  dislocation  density  at  the  surface  which 
means  the  existence  of  a  hard  surface,  but  if  the  sample  is 
left  for  a  long  enough  period  for  recovery  to  take  place 
such  that  a  reversed  dislocation  density  gradient  develops, 
then  a  hard  core  exists. 

Model  2 

The  two  dimensional  model  is  based  on  the  dynamic  mo¬ 
tion  of  parallel  edge  dislocations  moving  on  two  intersect¬ 
ing  slip  planes  at  an  angle  a,  as  in  Fig.  4.  For  FCC  metals, 
{111}  slip  planes  intersect  at  a  ■  70°.  For  each  dislocation 
present,  the  interaction  with  all  other  dislocations  and  their 
images  on  both  planes  have  to  be  consiaered  (Fig.  5) .  The 
applied  stress  component  on  the  secondary  slip  plane  is 
determined  from  the  Schmid  factor  which  corresponds  to  a 
maximum  probability  of  secondary  slip.  A  number  of  static 
interaction  for  various  dislocation  configurations  has  been 
studied.  The  dislocations  on  the  primary  and  the  secondary 
systems  seem  to  be  moving  in  a  cooperative  mode.  As  the 
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leading  dislocations  approach  the  junction  on  their  res¬ 
pective  plane,  they  encounter  increasingly  large  repulsive 
forces.  Therefore,  it  seems  that  if  both  systems  are  to  be 
operated  simultaneously,  they  will  block  each  other  at  the 
junction. 

If  the  inverse  pile-up  already  existed  on  the  primary 
plane  when  the  secondary  system  is  activated,  the  result  in¬ 
dicate  that  there  are  open-gate  and  close-gate  configura¬ 
tions  for  the  primary  dislocations.  In  other  words,  the  second¬ 
ary  dislocations  can  easily  move  through  the  junction  if  the 
primary  dislocations  are  preferentially  distributed;  whereas  in 
other  cases,  it  will  be  extremely  difficult  to  pass  it. 

Model  3 

In  FCC  metals,  edge  dislocations  on  the  secondary  and 
primary  slip  planes  are  never  parallel  to  each  other.  In 
order  to  represent  the  real  orientation,  a  three  dimensional 
model  is  needed  (Fig.  6) .  The  interaction  force  between  the 
primary  and  the  secondary  dislocations  is  thus  the  general 
case  of  a  non-parallel,  non-coplanar  dislocation  pair.  The 
correspondence  of  this  three-dimensional  model  and  the 
previous  two-dimensional  model  can  be  understood  by  looking 
at  the  viewing  plane  perpendicular  to  the  intersection  line 
MN,  (Fig.  7)  with  the  dislocation  penetration  spots  P  and  S 
considered  as  the  equivalent  two-dimensional  dislocation  lo¬ 
cations.  Therefore,  w.ien  the  two  dislocations  are  equi¬ 
distant  to  the  junction  point  0  and  are  on  the  same  side 
from  0,  they  are  actually  intersecting  each  other  along  MN 


somewhere.  As  the  dislocations  intersect  each  other,  an 

attractive  junction  should  have  formed  which  would  require 

2 

approximately  10  eV  to  break  them  apart  in  either  direction. 

An  interesting  fact  is  that  the  total  interaction  force 
between  a  general  dislocation  pair  is  independent  of  the  near¬ 
est  distance  between  them.  As  the  dislocations  approach  each 
other,  the  interaction  region  simply  becomes  more  localized 
with  higher  force  per  unit  dislocation  length  near  the  points 
of  smallest  distance.  An  attempt  is  made  to  simplify  the 
model  by  considering  only  the  length  of  dislocation  with  an 
interaction  force  per  b  greater  than  some  threshold  value. 
Since  the  total  force  is  a  constant,  as  the  dislocations 
approach  each  other,  the  effective  length  become  smaller  and 
the  effective  interaction  force  per  b  will  thus  increase. 

It  turns  out  that  the  effective  length  is  a  linear  function 
of  the  nearest  distance.  The  attractive  force  toward  the 
junction  on  both  planes  seems  to  agree  with  the  attractive 
junction  formation  arguments. 

The  plan  is  to  incorporate  these  three  dimensional 
results  into  the  two  dimensional  program.  For,  from  a 
computer  time  consideration  it  would  be  impossible  to  do 
dynamic  simulations  with  a  three  dimensional  model. 


TABLE  1 


AGo  (eV) 


Generation 
time  (sec) 


0.75  1.96 

0.85  27.3 
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Recovery 
time  (sec) 

2 . 8xl03  (1.3  hr.) 
3.0xl05  (4  days) 

3 . 9xl08  (10  years) 


1.05 
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Inverse  pile-up  of  dislocations  from  an  operating  source. 
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Fig.  3.  Dislocation  density  as  a  function  of  depth  beneath 
surface  before  and  after  recovery. 


.  Interaction  between  dislocations  on  the  primary 
and  secondary  slip  planes. 


